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Abstract 

The van der Waals forces across a very thin hquid layer (nanofilm) in contact with a 
plane solid wall make the liquid nonhomogeneous. The dynamics of such flat liquid 
nanofilms is studied in isothermal case. 

The Navier-Stokes equations are unable to describe fluid motions in very thin 
films. The notion of surface free energy of a sharp interface separating gas and 
liquid layer is disqualified. The concept of disjoining pressure replaces the model of 
surface energy. In the nanofilm a supplementary free energy must be considered as 
a functional of the density. 

The equation of fluid motions along the nanofilm is obtained through the Hamil- 
ton variational principle by adding, to the conservative forces, the forces of viscosity 
in lubrication approximation. The evolution equation of the film thickness is deduced 
and takes into account the variation of the disjoining pressure along the layer. 

Key words: Very thin films, interface motions, lubrication approximation. 
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1 Introduction 



The theory of thin liquid layers of a microscopic thickness is well understood 
(for a circumstantial bibliography, we may refer to the review article by Oron 
et al. [1]), but the motions of very thin liquid films wetting solid substrates 
are always object of many debates. In fact several problems appear: the liquid 
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in very thin layers is no more incompressible and the equation of motion is 
no more Navier-Stokes'; the concept of superficial energy related to a singular 
surface between gas and liquid layer has no more sense. 

Liquids in contact with solids are submitted to intermolecular forces inferring 
density gradients at the walls, making liquids strongly heterogeneous [23 . Of- 
ten, the fluid inhomogeneity in liquid-vapor interfaces was taken into account 
by considering a volume energy depending on space density derivative 
However, the van der Waals square-gradient functional is unable to model re- 
pulsive force contributions and misses the dominant damped oscillatory pack- 
ing structure of liquid interlayers near a substrate wall [1] . Furthermore, the 
decay lengths are correct only close to the liquid-vapor critical point where 
the damped oscillatory structure is subdominant [5]. Recently, in mean field 
theory, weighted density-functional has been used to explicitly demonstrate 
the dominance of this structural contribution in van der Waals thin films and 
to take into account long-wavelength capillary-wave fluctuations as in papers 
that renormalize the square-gradient functional to include capillary wave fluc- 
tuations [6j. In contrast, fluctuations strongly damp oscillatory structure and 
it is for this reason that van der Waals' original prediction of a density profile 
in hyperbolic tangent form is so close to simulations and experiments [7] . Con- 
sequently, depending of liquids and solids, a great number of different energy 
functionals may be proposed to model a liquid nanofilm in contact with a wall. 
To compensate the disadvantage of a special functional density to represent a 
thin film of liquid, we consider the most general case of any non-local density 
free energy functional and deduce a corresponding generalized chemical poten- 
tial. Then, we propose an equation of isothermal motions of flat nanofilms in 
contact with a plane solid wall. The classical chemical potential is obtained 
from the generalized chemical potential in the limit of homogeneous density. 
From the classical chemical potential it is possible to deduce the disjoining 
pressure value. The disjoining pressure exists only when the liquid is strongly 
nonhomogeneous [SfS] . The thin film is driven along the substrate by the dis- 
joining pressure gradient depending on the the layer thickness. 
We must emphasis that the derivation of governing equations is different from 
other approaches we can found in the literature. For example, in Ref. [Tj, the 
equation of liquid film dynamics is derived in the case of incompressible liq- 
uids. The disjoining pressure is ad hoc added as a formal additional pressure 
contribution. 



2 The disjoining pressure for horizontal liquid films 

The following experiment explaining the disjoining pressure concept is care- 
fully described in [9j. The liquid bulk of density was submitted to the 
pressure Ph by means of a very small bubble of radius R attached to a solid 
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wall (Fig. 1). 




Fig. 1. The bubble method of determining the disjoining pressure isotherms of wetting 
films (From [9], page 330). 

A thin film of thickness h separates a flat part of the bubble surface and the 
solid. The vapor bulk over the layer has a density p^^ and a pressure p^^- The 
difference between the two bulk pressures is 



where a is the surface tension coefficient of the bubble; the curve Il{h) obtained 
by changing the bubble radius is the so-called disjoining pressure isotherm. 

Without repeating the main results of Ref. [9] related to thin liquid films, we 
enumerate only the properties we use to describe thin films in contact with a 
solid wall. 

We consider the physical system at a given temperature 9 and suppose that 
the film is thin enough such that the gravity is neglected across the interlayer. 
The hydrostatic pressure in a thin liquid interlayer included between a solid 
wall and a vapor bulk differs from the pressure in the contiguous liquid phase. 
The forces arising during the thinning of a film of uniform thickness h produce 
the pressure Il{h) which is equal to the difference between the pressure p^^ on 
the interfacial surface which is - following the expression given by Derjaguin 
- the pressure of the vapor mother bulk of density p^^, and the pressure pb in 
the liquid mother bulk of density p^ from which the interlayer extends 



At equilibrium, this additional interlayer pressure n(/i) is called the disjoining 
pressure. 

The conditions of stability of a thin interlayer essentially depend on phases 
between which the film is sandwiched. In case of a single film in equilibrium 
with the vapor and a solid substrate, the stability condition is expressed as a 
form proposed in Refs. [lO], [Hj Chapter 4 : 



U{h) 



R 



n(/i) = Pv, - Pb- 



(1) 



dU{h) 
dh 



< 0. 
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Let us also notice an important property of a mixture of fluid and perfect gas: 
the total mixture pressure is the sum of the partial pressures of components 
and, at equilibrium, the partial pressure of the perfect gas is constant through 
the liquid- vapor-gas film where the perfect gas is dissolved in the liquid phase. 
Calculations and results are identical to those presented in the following sec- 
tions: the disjoining pressure of the mixture is the same as for the fluid without 
the perfect gas when a thin liquid film separates liquid and vapor phases (see 
Appendix). 



3 Equation of motion of thin liquid films 



In thin liquid films the density is strongly inhomogeneous. At temperature 6 
and for a given value /Xg of the chemical potential, the free energy of a liquid 
film in contact with a solid is a functional of the density p [12] . The free energy 
is written in the form: 



F = J j j^i'^ -pfiQ)dv + J J^(j)ds. 



In this expression D denotes the volume of the film and S the boundary of 
D. The boundary may be a free surface (as a liquid- vapor interface bordering 
the film) or a solid-fluid surface. The gravity effects are neglected. Moreover, 

a = pe = K [p] (3) 

denotes the free volume energy of the fluid and 

= L [p] (4) 

denotes the surface energy. Energies a and are functionals of the fluid den- 
sity p (for example a and are both functions of the density and successive 
gradients of the density). 

The equation of equihbrium is obtained by using the Gateaux derivative of 

r: 

n — /io = or grad = 0, (5) 
where fl is the variational derivative of a, 

We call Q generalized chemical potential of the fluid. For examples, the fol- 
lowing expressions of the volume free energy may be found in the literature: 

• The grand potential in molecular dynamics [5] 

ll I X^'^(l^-yl)pWp(y)^v (7) 



a 
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Then, 

The free energy for a generahzed van der Waals theory of capillarity [SlUSfP 



q; = q;(p, grad p, . . . , grad" ^p,0). (8) 

Then, 

where grad^ and divp denotes the gradient and divergence operators iterated 
p times. 

The expression of the surface free energy functional is expressed as an expan- 
sion in Ref. [Uj: 

, , 72 2 dp ^ dp . 

where 7^ = 1, ... ,4), are positive constants associated with intermolecular 
potentials between fluid and solid, the fluid density p is calculated at the wall 
and d/dn denotes the normal derivative to the wall. 

The solution of the equilibrium problem ([5]) with boundary conditions coming 
from the variation of functional ([21) determines an explicit form of the disjoin- 
ing pressure. An explicit example of such a calculation based on the van der 
Waals square-gradient functional 

a = Wip) + ^\Vp\\ 

where W is the volume free energy and A is a constant, is given in Ref. |15j . 

In the general case of conservative motions, the governing equations are ob- 
tained through the Hamilton principle of stationary action [TBIT7] : 
Let (t, x) be the Eulerian coordinates of particles. A particle is labeled by 
its position X in a reference space Dq. Let D{t) be a corresponding material 
volume. The motion of a continuum is a diffeomorphism from Dq onto D{t): 

x = v2(t,X), X = x(t,x), (pox = id. 

The equation of conservation of mass is 

pdetF = p„(X). F=?^. (9) 



5 



where pg is the density defined on Do, and F is the deformation gradient. Let 
us consider a one-parameter family of virtual motions 

X=$(t,x,r/) with $(t,x,0) = x(i,x) 

where r/ is a scalar defined in the vicinity of zero. The Lagmngian displacement 
associated with this family is P^16|18] 

5X=-(t,x,0). 

The displacement vector (5X is naturally defined in Eulerian coordinates. The 
Hamilton action is 



to 



Adv\ dt 



D{t) 



to 



>dS\dt 



where 



is the Lagrangian and V is the velocity vector . 

We consider virtual displacements (t,x) G [to,ti] x D(t) 

a neighborhood of S(t). From Eq. (Q we get : 



5X vanishing in 



6p = p divo 5X 



detF 9X 



6X 



where divg denotes the divergence in Dq. The definition of the Lagrange co- 
ordinates X implies 



9X(t,x)^ 9X(t,x) 



Then, 



(9x 



dt 



0. 



d6X^^ 9X,^, d5X ^ 
V + — 5V + — - =0, 
ax ax at 



or, 



6V 



d6X 
"dT 



where 



- = - + Vgrad 

denotes the material time derivative and the transposition. If we denote by 

= - - we obtain 
2 



5a 



dv 



to 



Do 



Po 



dt 



gradj N > 5X dvo 



dt 



dt 
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where gradg is the gradient in Dq. 
• The Hamilton principle of least action states 
V (t,x) G [to,ti] X D{t) — > (5X, (vanishing in the neighborghood of S(t)), 

Sa = 5 ri f f f Adv] dt = 0. 



Ito \J J JD{t) 

Consequently, the equation of motion is 

d (v^f) 



dt 

Let us note that 



grado A^. 



dt \ dxj 

where F is the acceleration vector. For a conservative isothermal motion, the 
equation of motion is 

r + gradr] = 0. (10) 

In nonconservative cases, to take into account dissipative effects, we simply 
introduce a symmetric viscous stress tensor cr^ in the second member of Eq. 
( fTOl) : the equation of motion becomes: 

^ + grad^] = -{div o-J. (11) 
P 

We consider a viscous stress tensor in the form 

cr„ = Ki{ tr D)1 + 2k2 D 

where 

D = -(— (^Y^ 

is the rate of deformation tensor and ni, K2 are the viscosity coefficients. 
Let us notice the important property: the gradient of the generalized chemical 
potential Q of thin layer is the only driving force of the liquid. 

Let us remark that Refs. [1.8J considered the liquid motion governed by the 
Navier-Stokes equations. The thermodynamic pressure was simply replaced 
with the disjoining pressure and the sign of the pressure gradient was changed 
in its opposite. An adherence condition was prescribed at the free liquid- vapor 
boundary layer. 
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4 Motions along a liquid nanolayer 



We consider a horizontal plane liquid interlayer contiguous to its vapor bulk 
and in contact with a plane solid wall. We use an orthogonal system of coor- 
dinates such that the z-axis is perpendicular to the solid surface. The liquid 
film thickness is denoted by h. 

When the liquid layer thickness h is small with respect to longitudinal dimen- 

f h \ 

sions L along the wall ( ^ ~ ^ ■'^ 1 ^^'^ inertia effects are negligeable, it 

is possible to simplify Eqs ffTTl) (approximation of lubrication [I9]). We denote 
the velocity by V = {u,v,w)^ where {u,v) are the tangential components. 
Due to the conservation of mass, du/dx, dv/dy and dw/dz are of the same 
order. The main terms associated with second derivatives of liquid velocity 
components correspond to d'^u/dz^ and d'^v / dz^. 

In the viscous stress tensor cr^ = Ki(divV) 1 + 2 K2 -D, the order of magnitude 
of partial derivatives implies 

CTt, ~ 2 /t2 D. 

As in ^20j, we assume that the kinematic viscosity coefficient u = 1^2/ P depends 
only on the temperature. Consequently, 

(1/p) div cr„ ~ 2z/ { div D + D grad [Ln(p)] } . 

In the liquid nanolayer, the hquid density variation \{p — Ph)/ Pb\ 1 and 
consequently D grad{Ln (p)} is negligible with respect to div D. Then, 

(1/p) div cr^ ~ 2 z/ div D. 

Hence, in the lubrication approximation the liquid nanolayer motion verifies 



grad f2 = z/ AV with AV 



(12) 



We denote by /i(p) the phase chemical potential obtained from Q for homoge- 
neous densities p; and of a plane liquid-vapor interface without solid wall 
and having the same value in the liquid and vapor bulks (e.g. p (p;) = p(p^) = 
0). In presence of a solid wall, the generalized chemical potential fl takes a 
bulk value p{ph) = p{PvJ = A^o (see Eq. ([5])). The value p^^ is the value of the 
vapor density at z — > +00, and p^ is a value of the mother homogeneous fiuid 
bulk in equilibrium with the vapor [9] . The z-component of Eq. (fT2l) yields 

J'^^ = and consequently Q[p] = p{pi,). 

To each value of p^ (different from the liquid bulk density value p; of the plane 
interface at equilibrium) is associated a liquid nanolayer thickness h. Then, 
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the functional fl[p] is a function of h: Q[p] = jiih) with = n^p^^ih)). The 
tangential components of Eq. (fT2|) yield for the velocity u ={u^vY' : 

(9^ u dp 0^ u 

grad/i(Pb) = ^ Q^^''^'^Pb= ^Q^^ (13) 

where grad denotes the two-dimensional gradient with respect to {x,y). 

Let us note that a liquid can slip on a solid wall only at a molecular level [21] . 
The longitudinal sizes of solid walls are several orders of magnitude higher than 
slipping distances and consequently, in a macroscopic mechanical model along 
the wall, the kinematic condition at solid walls is the adherence condition: 

uL=o= 0. 

From the continuity of fluid tangential stresses through the liquid-vapor inter- 
face of molecular size bordering the liquid nanolayer and assuming that vapor 
viscosity stresses are negligible, we get 



= 0. 

z=h 



du 

dz 



J^u= i^^z'^ -hzj giadpip^) . 
The mean spatial velocity u of the liquid in the nanolayer is 

u = — / u az. 



Consequently, Eq. (fT3|l implies 

-1 



h Jo 

and consequently, 

_ h'^ 
z/ u = -y grad/i(pb)- 

The chemical potential is homogeneous in the bulk; then, 

, / N dp jPb) 

dp[pb) = , 

Pb 

where p (p^) is the pressure corresponding to the bulk density p^. Hence, 

dp dp^ 1 9p(pb) 9pb , , 

grad ^(P.)-J^^^ g-d h ^ - — — grad h . 

The pressure p^^ = P (Ps^b) t^ie vapor bulk is constant and the disjoining 
pressure is 

U{h) =p(p^^ -p{Pb)- 
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Then, 

grad /i(pb) = gradn(h) , 

Pb 

and we obtain 

Xfe u = y grad n(/i), (14) 

where Xb = Pb^ denotes the dynamic viscosity of the hquid. The mean hquid 
velocity is driven by the variation of the disjoining pressure along the solid 
wall and the film square thickness. 

Eq. iHM is completely different from the classical hydrodynamics of films: 

Indeed, for a thin liquid film, the Darcy law is u = ~K{h) gradp, where p 
is the liquid pressure and K{h) is the permeability coefficient. In Eq. f|T4|) . 
the sign is opposite and the liquid pressure is replaced with the disjoining 
pressure. The mass equation averaged over the liquid depth is : 



|(/%rf.)+div (/%udz)=0. 



Since the variation of density is small, the equation for the free surface is 

dh , ^ , N 

— + /i div u = 0. (15) 

Replacing (IT^ into (ITSl) we finally get 

dh ( \ 

— + divh— gradn(h) =0. (16) 
ot \^Xb J 

Eq. fll6p is a non-linear parabolic equation with a good sign, it ^<o. 

oh 

This result is in accordance with a physical stability argumentation about thin 
hquid layers obtained by [ 91ITT] . 



5 Conclusion 



In this paper we use the variational principle of stationary action to obtain 
the equation of motion of a liquid in a very thin plane layer. The result is 
obtained whatever is the functional of free energy in the liquid and the dis- 
joining pressure gradient is the driving force of the film motion. This result 
is obtained in a direct way without any phenomenological assumption as it 
is done in the literature where the disjoining pressure of a thin film is added 
formally in the equation of motion of incompressible liquids. The disjoining 
pressure expresses the excess free energy of the layer with respect with the 
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superficial energy of a fluid interface (superficial tension) [22] . Consequently, 
in the equation of motion, the superficial tension has no raison to be taken 
into account. This result can be different if some curvature appears in the film 
as in the case of a micro droplet. We have also obtained the evolution equation 
for very thin isothermal flat liquid films. The thin film is driven along the sub- 
strate by the disjoining pressure gradient associated with the layer thickness. 
The equation of the film thickness evolution is a nonlinear parabolic equation 
having a good sign when the thin liquid layer is stable, the instability being 
probably associated with the rupture of the film along the substrate. 
In appendix, we will see that the results are unchanged when the layer is 
constituted with a mixture of fluid and gas. 
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6 Appendix: Extension of results for a fluid mixture 

Let us consider a mixture of the perfect gas and the non-homogeneous hquid. 
A mixture of the perfect gas and the saturated vapor of the liquid rises above 
the thin layer. To model the mixture at temperature 6, we consider a volume 
free energy in the form [23] : 

pe = a + P, 

where a = K[pj-] is defined as in Eq. and P = H{pg) is the free volume 
energy of the perfect gas. The mixture density is p = pj + Pg where is the 
fluid density and Pg is the perfect gas density. 
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In conservative case, without diffusion between constituents and neglecting 
the gravity forces, the equations of motion are [23] : 

^;+grad^][p^] = 0, ^^^^ 
Tg + grad u{pg) = 0, 

where F, {i = /, g) are the accelerations of components, ^[Pf] is the generalized 
chemical potential of the fluid (Eq. (1141) ) and uj{pg) is chemical potential of 
the perfect gas. 

At equilibrium, the thin layer is invariant along the solid wall and Eq. ffTTl ) 
yields Pg = PgQ where PgQ is a constant in the mixture. The total pressure in the 
mixture bulks is the sum of the partial pressures of constituents: p = Pf + Pg- 
The partial pressure of the perfect gas is constant in the domain occupied by 
the mixture and the disjoining pressure is: 

n = {Pfv - Pg) - (Pfb - Pa) = Pfv - Pfb- 

Hence, the disjoining pressure can be reduced to the disjoining pressure for 
a single inhomogeneous fluid. If we take into account the liquid viscosity and 
neglect the perfect gas viscosity, Eq. (fT71i ) becomes: 

Tj+grad n[pj] = ujAV, 

where z^/ is the fluid kinematic viscosity. This result can be extended to multi- 
component cases. 
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